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INTRODUCTION 


It is universally assumed in the text-books on 
suspension bridges that the shape of the curve 
assumed by the cable isa parabola on the suppo- 
sition that the weight of the structure and the loads 
are uniformly distributed in a horizontal direction. 
Neglecting for a moment the moving loads, the 
weight of the structure is made up of the weight of 
the cable and the suspenders, and the weight of the 
roadway and the stiffening strusses, it being assumed 
that the stiffening is not of the braced-chain type. 
We could simply classify this as the weight of the 
cable and the weight of the roadway, the former 
being uniformly distributed along its length and the 
latter horizontally. It is well-known in statics that 
if the former be inappreciable and can be neglected 
the curve of equilibrium is a parabola; while if 
the latter be neglected in comparison with the 

cable weight the curve is a catenary. It is hardly 
necessary to point out that the curve will still 
be a parabola even if the former be not neglected 
but taken to act horizontally. 1 propose to examine 
in this paper how far this assumption Is justifiable 


* Paper presented ‘at the 26th Annual Session of the Mysore 
Engineers Association held in March 1934. 
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and what is the error made by the assumption in 
calculations relating to suspension bridges. In view 
of the fact that every book, either on applied 
mechanics or suspension bridges, makes this 
assumption without any attempt at justification, it 
is worthwhile to investigate how far it is correct. 
Before proceeding to this discussion, it might be 
observed that greater errors are likely to result from 
this assumption, the greater the weight of the cable. 
The following *table giving a comparative statement 








Net weight of Gross fixed 
Name. Span in ft. chain per ft. load per ft. 
run. of span. 
Clifton AUg.23 1560 3171 
Freiburg 870 167 760 
Niagara O4les 820 2032 














for three suspension bridges shows that the cable 
weight is quite appreciable when compared with the 
gross fixed load and it may consequently be of some 
practical use to undertake this investigation. 
Further, it is also evident that large values of the 
ratio of sag to span introduce greater errors in the 
assumption and in practice this ratio varies from 
0.2 to 0.05 such that it might be as high as 1/5. 


§ 1. THE GENERAL DIFFERENTIAL EQUATION 


To determine the form of the cable for any 
loading, consider vertical loads applied on a cable 
suspended between two points. The cable will 
assume a definite polygonal form determined by the 

“ Table from Rankine : Civil Engineering. p. 582, 4 
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relations between the loads. The end components will 
be inclined and will have horizontal components H 
Simple considerations of static equilibrium show 
that H will be the same for both end reactions and 
will also equal the horizontal component of the 
tension in the cable at any point. H is called the 
horizontal tension of the cable. It is easily shown 
that the cable curve is the bending moment diagram 
for the applied loads. If the loads are continuously 
distributed the funicular polygon becomes a conti- 
nuous curve and if w be the load per horizontal 
linear unit at any point having the abscissa x, we 
can obtain the equation for the vertical shear at any 
section x of the span and on differentiating this 
equation we,obtain, 
Hoy 
dx” 





— Se ed. 


and on taking the y-axis vertically upwards this 
becomes 





ry aNao 
Hance 


which might be taken as the differential equation* of 
the equilibrium curve in the general case. Let g, be 
the weight of the cable per unit .length (supposed 
uniform) and w. the total uniform moving and dead 
loads per, horizontal foot ; equation (1) then becomes 


da’ 
Ho = go secX + Wo (2) 


where « is the inclination of the cable to the horizon- 
tal at any point. Eqn. (2) is the general differential 
equation of the problem. If go be assumed to be 








*“ See for e.g. W. H. BURR: Suspension Bridges p. 61. 
Eqn. (13) or JOHNSON, BRYAN and TURNEAURE: Modern 
Framed Structures. Part II, pp. 187-88. We shall use 
throughout the notation of the former book, 
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uniform horizontally instead of along the curve, we 
can reduce the equation to 


2 
H a= Ww (w=Wot go) 


and taking the initial conditions at the lowest point 
of the cable to be: 


x =0, =0, y = O, this equation gives 
x 
2 2 
Hy: eae a Ae 
y Ory 2H 
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which is the standard equation to the parabolic cable. 
The tension T at any point is given by 


T =H Sec « = HA, * (22) 
a 


If the points of support are at the elevation y =f 
and the span = /7, H takes the value 


wt we / A 
H="j and T = iV +(¥) 
If L be the entire length of the cable we have 


—9 [ I (22 214, 1 ; . 
e =2 [ Ir ee) | dye (2u+sin h2u) * 
c= L/ 64f2 2 $ 

2 fi [ 1 oe ] Ke 


Developing the radical by the binomial] theorem to 
three terms and performing the integration indicated 


L=1| 1+ 8 F432 i 
= a Saad, 5 
The equations for y, H, T, L are the fundamental 
equations for the parabolic cable. 











“n = f/land u is defined by sin hu = 4 n.%; This elegant 


form for L is to be found in Steinman; Suspension Bridges p. 5 
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If in equation (2) we neglect w, we obtain 


2 
HS = gy sec & =g.A/ 1+ (22) 
This is the differential equation of a catenary and 
its solution gives 


> Oats | Cx —Cx 
ee (. ap ae) where c= g./H 
and we can derive corresponding equations for H, T 
and L as in the parabolic cable. 


I now consider the general differential equation 
(2) where neither w, nor g. can be neglected (it 
must be observed that assuming a uniform horizon- 
tal distribution of g, is mathematically equivalent to 
neglecting it) and obtain the exact mathematical 
solution. In the next article, I derive an approximate 
form of the solution suitable for practical use and 
enabling us to compare the formule obtained for a 
strictly parabolic cable. 


We shall write (2) in the form 


HE =¢. sec © + wo (replacing « by © ) 


ae 
aot =8 p= 
or Caan Be Oot bi (« We =) 


and assume, as always realised in practice, that b>a. 
Putting p=dy/dx, this further reduces to the first 
order differential equation: 

dp _ 1+p2+b. 3 

dg_ / (3) 

Since p = tan ©, dp = sec? ©. dO and (3) can be 
written as 


_sec?6.d6_)7,. d © Oa 
b + asecO ; cos © (a+b cos QO) 





dQ ibd 





2 =adx 
eye ay e eeey operat + 
and integrating we have 
ax = log (sec 8 + tan O)- 
b Pap b?—a? sin© 
Phe te 6 (6 +a cos a/ a (4) 


Nias log | aoe 


the constant of integration being given by x = 0 for 6 = 0. 


a+bcos O 





Again oP a 1 + p? + b can be written as 
2 





__pdp ce tan © sect Od |, 
afi +p? + d a asccO+5 —% 


This can be easily integrated and the result is 


+ 0d 
eee Hi} 4 735] (5) 
a®* y = a(sec 0-1) b log Stag ears 


taking y = 0 for 6 = 0. 


Equations (4) and (5) give the parametric equations 
of the curve of the cabie in terms of ©, but they are 
in quite an unmanageable form for any practical 
purposes. Before deriving approximate equations 
from these, it is instructive to see how these reduce 
to the parabola and catenary for thecases a=O and 


b=0 respectively. 

Case b= 0: Equation (4) takes the simple form 
ax = log (sec © + tan ©) | 

and equation (5) reduces to 


ay = sec @ —1, since the argument of the 
logarithm in equation (5) can never become infinite. 
Eliminating © between these equations we have 
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= ax —ax ) 
ae a 


which is the equation to the catenary 


Case a = 0: In this case equations (4) and (5) reduce to 
identities but by dividing them throughout by a and 
a’ respectively and applying de L’Hospital’s Rule 
for the evaluation of indeterminate forms, they 
take the form 


= tan © 
1 
and y =p tan? © } 


: bx? 
which give y ae , the equation to the parabola. 


§ 2. APPROXIMATE SOLUTION OF THE GENERAL 
DIFFERENTIAL EQUATION. 


The differential equation (3) is equivalent to the 
two equations 


dy _ 2 


dx Wied 


ee 2 EE eee 
dx 





where zg = sin ©. Assuming for yand z two power 

series * with undetermined coefficients, using the 

initial conditions x=0, y=0, z= O and determining 

these coefficients so as to satisfy the above two 

equations, we deduce for y the expression 

_(atb) 2 | abex' 5 
eee Sag. it ' (5) 

“This method is suggested for approximate solutions of 


differential equations in RUNGE UND KONIG: Numeris- 
ches Rechnen. § 87. pp 321—23 
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As a first approximation we can retain only the second 
term in (5), and the equation then illustrates how 
the assumption of non-uniform horizontal distribu- 
tion of g. affects the shape of the parabolic cable. 


I shall now indicate an alternative method of 
deriving equation (5) which is more illuminating. 
Equation (2), vez. 


d 
ob = aN 1+ pt b (2) 


gives us the parabolic cable when we write @ in 
place of a v1 + p? which is equivalent to putting 
p= 0 2e., treating the cable as straight. We can 
therefore obtain a better approximation by putting 





for p Se) its value as obtained for a parabolic 


cable. From § 1, the equation 


b x2 
y= aye gives p = =bx and hence (2) can be 


approximated to 





dp _ 22) 
ae a(l+obx)t+ 5 

edits | ab Nea | 

1.€., ——— = @ 1+ aaa eee a eC 
dx* fA 8 hae 


on expanding by the binomial theorem, since bx 
giving the slope at any point can be assumed, for all 
practical purposes, numerically less than unity, and 


the expansion being therefore valid. Retaining 
only two terms, this gives 

d2y b2 x2 

Aa ae ae ee i 


Integrating once and using oY =0 iss 7) 
x 
dy 


—~=(a + b) x + ab? xt 
dx 6 
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Another integration, using y = 0 if « = 0 gives 
(a + b) x2 ab’ x‘ 
a ee eee 5 
2 24 6) 


* 





which is the same * as equation (5). For use in 
practical calculations it is preferable to express this 
equation in terms of fand/; but we have only one 
condition expressing that (5) passes through the 
points of support ( + 7/:, f)and_ this does not enable 
the determination of both aand 6b in terms of these 
quantities. Let us, however, write a—=6n, so that n 
is the fraction of the cable weight to the horizontal 
load; equation (5) then reduces to 
_ b(n+1) x? nb! x 
2 24 


Putting x=1/2 & y=f, we have 
=o (n+1)?? 4 nb! e 
8 384 


bl? _nbl? eS 
Ba f— eo = (10 +i rs 


sf 


This equation gives a cubic in 6 from which we can 
find b in terms of n, f and 7 which might be supposed 
to be known quantities. We can solve this cubic 
equation approximately. 


2 
Writing 2 =, eqn. (6) reduces to 





Bate ( 492 \- A 
f no€ 1+ 3 2 , a cubic ind 


For small values of , a first approximation is 


Ft aor 
ta PD ales 





*~ See W. H. Burr: ibid p. 69 where a similar method is 
adopted in a different connection, 
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which is a well-known expression in a parabolic cable. 
To obtain a second approximation let us substitute 
this value ) = fin the right hand side and obtain 





4 =f—nf (1 + iF )=s (1-n— SF 





; Tah _.. 4nf? 
1.6, 6 = 7 1l-n it) (7) 


8 8 
nb x 


(8) 





Alen: “=p (aut ee 
Oy 


6 being given from (7) in terms of n, f and 1. 


Let us now calculate the length of arc of the 
entire cable as modified by the above formule and 
compare the same with the original expression 


Laji+S HF 
3 12 Dea a 


obtained in § 1. Denoting the modified length as 
L’ we have 


’ l/s | dy\2)4 < 
iz =o f¥ il + (2 |? dx :@ being given from (8) 


/ 2 8 2 
= [D+ { (a+) et Oe } ]i« 


Expanding by the binomial theorem and retaining 
terms up to ,*‘, this reduces to 


v=o ffi 4 (a+b)? { (a:tBh tat so 8) ) aa 
: 2 24 j 


=r fis (a+ OV , (a+b) § 4ab2—3 (a +0)" } 
24 1920 
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on performing the integration indicated. Now, to 





substitute @ = bu and b= 4 eta be 


372 p in this ex- 


pression, we can write 


64 fF 2 4 nf2\2 
2 tI TS! 
b rz (: n ie) 


= 64 f? 8n ae : 
7! (: Zn ane , retaining only first power of n, 





Similarly '= aoa ye eee as 


Bh 


2 
Now L'=/ fin He gels (1 + 21) my [+n) 


(3 n°+3 n pene, 


= fi+ ass n)— ve 2 )} 


neglecting terms in 22, n°, etc. 








or path i a de ea wee ) 
| B45 9 1920)..12 so) | 


Substituting in this the values of 5? and 64 written 
above, we get finally 








~~ Bera acre) fag fF! BIZ ht 
oS aapeae 7B) hag tar tps 


~ Hence L’ — L=1! 8 n 4/14 (neglecting the last term) 
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In practice it is found that all ratios of #/: will 
be between .2 and .05 and the following table * 
gives values of L/l corresponding to different values 


of #/1. 








f[l L/l 
2 1.09643 
175 ) 1.07566 
15 | 1.05676 
125 | 1.040104 
1 | 1.026027 
075 1.014797 
05 1.006627 





ernest 
* Vide Burr : ibid. (p. 73) 


UHL 148, 7 

l ese 
gives us an indication of the error committed in 
treating go as a uniform horizontal load. Taking from 
the table in the Introduction the case of the Freiburg 
suspension bridge for which n= (approximately) let 
us calculate the deviation for the case f/l=.2; we 
easily find 

E = 1 ="005262 

such that the corrected value would be 1.10169, 
showing that although the modification required is 
not great, it is of the same order as the difference 
between the values of L/z calculated as in the above 
table and the value calculated from the formula t 
exact for a parabolic cable, viz. 


1 
{}J==— (9 .j 4) 
L/l én (2M + Sin h2u ) 





The above expression 











t Vide Burr: ibid, (p. 7) 
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§ 3. OTHER DIFFERENTIAL EQUATIONS 


In the previous investigation we have considered 
the weight of the cable alone and neglected that of 
the suspenders. Let g’ be the weight of the suspenders 
per square unit of vertical plane of the cable so that 
gy is weight per linear horizontal foot of that part of 
the suspenders above the vertex of the cable. The 
part of the suspenders below the vertex forms a 
uniform load over the span and might beassumed to 
have been included in W,. The general differential 
equation becomes 





2 
“= NOa tpt bDiog 





1.€., a Ee AAG + p? +b + cy (c=g'/H ) 


Making # and y the dependent and independent 
variables respectively, this can be reduced to the first 
order differential equation. 





pe =a n/ +p>+b+cy 


of which an exact solution is complicated. It can 

however be solved by approximate methods by 

substituting for y and 4 the expressions as given 
dx 


by the parabolic cable. 


2 
Putting y= and ie equation (10) reduces to 
2 dx xa 


tnt Ag A ioe dee 
oy =a A 1+ 6? x 2 


=a ( 1+2¥)+o+ a 


=(a+5)+ "= (b +0) 
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Successive integrations give 


PS (teh) apse ee eee 
dx 6 

To (gab) te een ee 
errs Log re 


which expresses the deviation from the parabolic 
form. 


I now consider, finally the case where the cable 
is not necessarily uniform. Let g, be the weight per 
ft. run at the vertex and g, at any point where 
inclination of the cable to the horizontal is 6, 
supposing that the normal section of the cable at 
any point is proportional to its total stress, we have 


& x=. SecO x SecO =g, Sec? ©. 
and the differential equation becomes 


ig Pe 

— =a Se?O6+b+cy (11) 
dx? 

Putting c= °F and y= +f x2 


this can be written 


San dee 6472, 4cf 
Tene Pn Saati x2 





2 
giving oe (a+b) +kx* 


Bie e ei 

h ee ee ie et Le 

where k a a +e) 

An exact solution of equation (11) is a matter of 
some difficulty. 


I hope to return in a further paper to more 
detailed discussions of equations (10) and (11) and 
also discuss, if any, the bearing of the results of § 2 
on the problems of stiffening of the structure and the 
deflection of the bridge system, 
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